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The general exact solution describing the dynamics of anisotropic elastic spheres supported 
only by tangential stresses is reduced to a quadrature using Ori's mass-area coordinates. 
This leads to the explicit construction of the root equation governing the nature of the 
central singularity. Using this equation, we formulate and motivate on physical grounds 
a conjecture on the nature of this singularity. The conjecture covers a large sector of the 
space of initial data; roughly speaking, it asserts that addition of a tangential stress cannot 
undress a covered dust singularity. The root equation also allows us to analyze the case 
of self-similar spacetimes and to get some insight on the role of stresses in deciding the 
nature of the singularities in this case. 
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I. Introduction. 

The present status of research in black hole formation and cosmic censorship is quite 
intriguing. Indeed, both analytical results in dust collapse (see e.g. Jhingan & Joshi 1997 
and references therein) and numerical results in scalar field collapse (see e.g. Gundlach 
1997 and references therein) indicate the existence of a critical behaviour governing the 
formation of black holes or naked singularities. In principle, this behaviour should be the 
"remnant" of some hypothesis of a - still unknown - cosmic censorship theorem and, as 
such, should be related to the properties of the collapsing matter like fulfilment of energy 
conditions and local stability. If we want to understand the physics of such phenomena 
in the case of "ordinary" matter (i.e. not boson stars), we are enforced to approach 
analytically the gravitational collapse of non pressureless matter, since the dust equation 
of state is "trivial" from the viewpoint of matter properties: all energy conditions, causality 
conditions and stability conditions reduce to positivity of energy. However, this problem 
is extremely difficult if approached in full generality (see e.g. Joshi 1996 and references 
therein). 

Recently (Magli 1997, to be referred towards as [I]), we discussed a class of solutions 
of the Einstein field equations describing spherically symmetric, non-static elastic spheres 
supported only by tangential stresses (elastic matter is of interest in strongly collapsed 
situations; for instance neutron stars typically have solid regions - see e.g. Haensel 1995). 
Previous investigations on systems having vanishing radial stresses trace back to Einstein 
(1939) and Florides (1974) in the static case, while non-static models have been considered 
by Datta (1970), Bondi (1971), and Herrera & Santos (1995) (for a general review on 
anisotropic systems in General Relativity see Herrera & Santos 1997). 

From the viewpoint of cosmic censorship, understanding the nature of singularities 
for such solutions can be considered as a first step toward a full understanding of spherical 
gravitational collapse with general stresses. Indeed the investigation of some particular 
models carried out recently by Singh and Witten (1997) already shows behaviours which 
can be drastically different from the dust ones. The purpose of this paper is to contribute 
to a program which should hopefully lead to a full understanding of the final fate of 
gravitational collapse with tangential stresses. 

The paper is organized as follows. The first problem we have to face with is the fact 
that, due to the presence of stresses, the comoving time differs from the proper time of 
the shells of particles. As a consequence, the first order "energy equation" which arises 
from mass conservation is coupled, for generic equations of state, to the equation giving 
the acceleration of the world lines of the particles. This coupling has the effect that 
the equation for null geodesies (whose behaviour near the central singularity governs the 
nature of the collapse) cannot be written in explicit form. Here, we overcome this difficulty 
(section II) using mass-area coordinates. These coordinates were originally introduced by 
Ori (1990) to obtain the general exact solution for charged dust. Here we obtain the 
general solution for gravitational collapse with non vanishing tangential stresses in a very 
simple form, in which only an integral remains to be performed. 

Using the line element in the new coordinates, it is possible (section III) to write 
the "root equation", namely the algebraic equation governing the nature of the central 
singularity, in explicit form. This equation depends only on the choice of the equation of 
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state and of the initial distributions of density and velocity. 

It turns out to be a quite difficult task to study this equation in full generality; it is, 
however, worth mentioning that a full understanding of the nature of the singularities for 
dust spacetimes has been achieved only very recently (Singh & Joshi 1996). To get some 
insights into the physical contents of this equation, we carry out a detailed examination of 
small deviations from the dust equation of state. This leads us to formulate a conjecture 
on the final fate of the collapse and to give a plausibility argument supporting it (section 
IV). Finally, in section V, we analyze the case of self-similar spacetimes putting in evidence 
some qualitative effects of the tangential stress on the nature of the central singularity. 

The paper ends with some concluding remarks in section VI. 



II. Mass— Area coordinates and general solution 

Recently (see [I]) we discussed a class of solutions of the Einstein field equations 
describing spherically symmetric, non-static elastic spheres supported only by tangential 
stresses. Using comoving coordinates, the line element reads 

ds 2 = -e 2v dt 2 + ^^-dr 2 + Y 2 (d9 2 + sin 2 ed<j> 2 ) , (1) 
where v and Y are function of r and t satisfying 

yV* = -1 + H + i + Z . ( 3 ) 

In the above formulae, F(r) and f(r) are the "conserved mass" and the "binding function" 
familiar from the Tolman-Bondi solutions, while h = h(r, Y) is the internal elastic energy 
per unit volume (it plays the role of equation of state of the material). For physical 
reasonability, F, 1 + / and h must be chosen as positive functions (positivity of mass, of 
g rr and of elastic energy, respectively). Moreover, the function h is severely constrained 
by the local stability of matter, which requires this function to have a minimum at Y = r 
(see [I] for details). If h is equal to one (more precisely, is a constant which may be 
rescaled to unity) there is no dependence on the strain: the material is a dust cloud and 
the line element reduces to the Tolman-Bondi one. This fact will be very important in 
what follows; we shall call dust limit of any equation depending on the choice of h the same 
equation written with h = 1. 

The energy-momentum tensor of the material is diagonal in the comoving frame and 
has only three non-vanishing components, namely the energy density e = — Tq and the 
tangential stress II = Tq =T?. These quantities are given by 
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n = He , (5) 
where the "generalized adiabatic index" H is defined as follows: 

In [I] the physical properties of these solutions are thoroughly discussed. We recall here 
that the metric is well behaved at the centre if the equation of state satisfies the "minimal 
stability requirement" (h has a minimum at Y = r) and the conditions 1^(0, t) = 0, 
/(0) = h 2 {0,0) - 1 hold. The behaviour of the energy density at r = is the same 
as that familiar from the Tolman-Bondi models, namely e is initially regular if F(r) is 
of the form r 3 F(r) with F(0) < +oo. Matching our solutions with the Schwarzschild 
vacuum is possible on any chosen boundary surface r = r&, provided that the value of F 
at Tb is identified with the Schwarzschild mass M. However, the transformation between 
comoving and Schwarzschild coordinates is highly non-trivial. The energy conditions lead 
to inequalities on the function H, and therefore to differential inequalities on the state 
function h. In particular wee holds if H > —1. Once this is satisfied, dec requires H < 1, 
while sec is satisfied if 7i > —1/2. Analyzing the behaviour of the function appearing at 
the right hand side of equation (3) for a fixed shell of particles (r = const.), one can give a 
qualitative analysis of the possible motions. In particular, it is shown in [I] that there exist 
physically valid models of oscillating elastic spheres as well as of finite-bouncing spheres. 

The metrics (1) are not completely explicit due to the coupling between eqs. (2) and 
(3); physically, this is simply a reflection of the fact that the comoving time differs from 
the proper time since the particles are not in geodesic motion. The two equations decouple 
only if h depends uniquely on7 (h = w(Y), say). In this case we have 



e 2v = 



v2 



1 

W 2 



1 


"1 + / 


(-¥)] 


w 2 


w 2 





and the function H defined in (6) reads 

w = 4^io gro( y). 

This particular class of solutions contains that discussed by Singh & Witten (1997), in 
which the tangential stress is proportional to the density. In this case the stress-strain 
relation has the "barotropic form" IT = ke with constant "adiabatic index" H = k, and 
the equation of state is w(Y) = Y~ 2k . 

A consequence of the above described coupling problem is that, in general, it is not 
possible to write explicitly the null geodesic equation in comoving coordinates. As will be 
recalled in the next section, this equation governs the nature (naked or black hole) of the 
central singularity, and it is therefore very difficult to investigate on censorship using the 
comoving frame. 



G. Magli, Gravitational collapse with non-vanishing tangential stresses II. 



4 



A system which is mathematically very similar to ours is charged, spherically symmet- 
ric dust (Vickers 1973). Indeed for such a system the mass is conserved and it is possible 
to analyze the dynamics in comoving coordinates; however the coupling between "times" 
does not allow explicit integration. To get rid of this problem Ori (1990) introduced a 
system of coordinates which removes the coupling and obtained the general exact solution 
for charged dust. Ori's system is obtained by replacing the comoving time t and the radial 
label r with the "area coordinate" R = Y(r,t) and the "mass coordinate" m = F(r); the 
mass coordinate, being conserved, is comoving (if denotes the velocity of matter, one 
has ViP 1 = uSfi where u := u R ). The line element in these coordinates has the form 



where A, B and C are functions of m and R. 

We are now going to show that Ori's technique can be applied also in the case of 
solutions with vanishing radial stresses. This is essentially due to the fact that the mass is 
conserved also in this case, and therefore gives an unambiguous "comoving label" for the 
shell of particles (for simplicity, the general solution is presented here in the non-charged 
case; however our results can be easily extended to the case of charged materials, as briefly 
reported in appendix). 

Since m is comoving, we have 



The energy density can now be written as 



where h is to be considered as a function of m and R, the quantity H is defined by 



and E is an arbitrary function of m corresponding to y/1 + f (we have introduced this 
notation in order to facilitate the comparison with Ori's 1990 paper). 



There are four (compatible) Einstein equations for the three unknowns A, B, C. We 
start considering equations G™ = 8nT™ G* = 8irT^ and G% = 8ttT%: 



ds 2 = -Adm 2 - IBdRdm - CdR 2 + R 2 (d6 2 + sin 2 9dip 2 ) , 



h 



(7) 



e = 



4ttuR 2 eVh 



H :=B 2 -AC , 






(9) 



(10) 
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Equation (8) can be integrated: 



A = H [ 1 - ^ 
R 



(11) 



Using B 2 = H + A/u 2 and equation (7), from (9) we get 



u = ± 



2m E 2 

1 + ir + ^ 



(12) 



and the metric function B can be written as 



B 



hu 



Therefore, we have solved for u (and thus for C) in terms of the arbitrary functions and 
expressed A and B in terms of a single unknown H. To complete the solution, we plug 
the above results in (10), obtaining 



h (\ 



,m 



,R E \U 



We have, therefore, reduced the problem to the calculation of an indefinite integral: 



y/H(m,R) = g(m) ± J G{m, R)dR 



(13) 



where the ± sign is the same as that of u, g(m) is an arbitrary function, and 

h 



G(m, R) := 



RE 



R (E 2 
1+ 2 (W 



,m 



_ 2m E 
+ 



(14) 



If h = 1 the above formulae give the Tolman-Bondi line element in mass-area coordinates 
(Ori 1990). 

It is easy to check that the remaining field equation G e = 8ttTq is identically satisfied 
once eqs. (8), (9) and (10) are. 



III. A laboratory for Cosmic Censorship 

In this section we use the general exact solution derived in the previous section to 
build up a "laboratory" for studying cosmic censorship. The key instrument which is 
needed in this laboratory is already known from the work by Dwivedi & Joshi (1994) and 
may be called "root equation" ; as recalled below, it is an algebraic equation arising from 
the behaviour of outgoing null geodesies near the singularity. In the present section we 
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construct this equation explicitly for the case at hand. We shall also derive the conditions 
for shell-crossings singularities in terms of an integral equation. 

3.1 Physical content of the arbitrary functions: initial data 

In order to approach the problem of singularities, we first need to identify the physical 
content of the arbitrary functions, so that regular initial data can be chosen. 

There are three arbitrary functions, namely the equation of state h = h(m, R) and 
the functions E(m),g{m). To understand the physical meaning of E and g observe that, 
physically, such functions must be related to the "initial distributions" of density and 
velocity (here quotation marks are due to the fact that we shall take care of the initial 
data always referring to the "original" - comoving - coordinates). Consider, therefore, 
regular initial data at some comoving time t (t = 0, say). We use the scaling freedom in 
the choice of the r coordinate to identify the lagrangian and the eulerian label initially, so 
that Y(r, 0) = r. In mass-area coordinates, to the equation Y(r, 0) = r corresponds some 
curve R = R (m), where R = F -1 and we are assuming the mass to be a monotonically 
increasing function. Introducing the initial distribution of velocity (V(m), say) from (12) 
it follows 

Tr2 , . E 2 (m) „ 2m 

The above formula gives the relationship between E and the initial velocity profile. 

The relationship between g and the initial data is, in general, quite complicated. To 
obtain it, observe that the following formula may be easily proved: 

The above equation evaluated "at t = 0" yields 



g(m) = ± 



R° m (m)h(m,R ) f R °^ 
V(m)E(m) J 



G(m, R)dR 



(16) 



(the ± sign is the same as that of u) . 

3.2 Shell-focussing singularities 

The energy density (4) becomes singular whenever Y(r, t) or Y'(r, t) vanish during 
the dynamics. Physically, such singularities correspond to those occurring in dust models: 
Y = corresponds to "crushing to zero size" (shell-focussing singularities) while Y' = 
corresponds to the shell crossing phenomenon: the world lines of the (shells of) particles 
intersect each other and the "lagrangian labelling" description breaks down. 

Contrary to what happens in the dust case, where shell-focussing collapse in unavoid- 
able, within our solutions there exist globally regular models of oscillating or bouncing 
back materials. However, equation (3) can be used (see [I]) to show that for any physically 
valid choice of the equation of state it is possible to choose initial data leading to continued 
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gravitational collapse and therefore to shell-focussing singularities (it is worth mentioning 
that the remark made by Singh & Witten (1997) that regularity conditions "explicitly 
disallow the formation of a singularity at r = 0" is incorrect: these conditions imply only 
regularity on the initial data surface and non-preferredness of the centre). 

To analyze the nature of shell-focussing singularities, we first observe that the relation 
Y(r,t) = defines a "singularity curve" t s (r); in general, different shells become singular 
at different times, and it is customary to call central singularity that occurring at r = 0. 
This singularity plays a distinguished role because it is possible to show that non-central 
singularities are always covered. For, notice that the shell labelled r becomes trapped 
at a time t t (r) such that Y(r, t 1 ) = 2F(r). For each fixed shell, consider the function 
Y(t) = Y(r,t). We have Y{f) = 2F(r) > Y(t s ) = 0, but dY(t)/dt is negative in a 
collapsing situation, so that Y(t) is decreasing and it must be t s (r) > t t (r). It follows 
that the shell becomes trapped before becoming singular, so that the singularity is covered 
(under certain conditions, it is possible to proof this result also in presence of non-vanishing 
radial stresses, see Cooperstock et al (1997) for details). 

The above argument does not work for the central singularity, at which F(0, t*) = 
2F(0) = 0. To study this singularity, we translate in the mass-area formalism and make 
use of the method developed by Dwivedi & Joshi (1994) and successfully applied to the 
dust case in a series of recent papers (see e.g. Singh & Joshi 1996, Jhingan et al. 1996). 

Consider the equation for radial, outgoing null geodesies in mass-area coordinates: 



dm C \ h 

This is an ordinary differential equation with a singular point at the central singularity 
R = 0, m = 0. This singularity is (at least locally) naked if there are geodesies starting 
at it with a definite value of the tangent. If no such geodesies exists, the singularity is 
not naked and (strong) cosmic censorship holds. To investigate the behaviour near the 
singular point, define 

R 



x 



2m c 



where a > 1/3. If the singularity is naked, there exist some a such that at least one finite 
positive value xq exists which solves the algebraic equation 

Xq := lim x = lim 



R,m^0 R,m^0 2m c 

Applying L'Hospital rule we have 



m x ~ a dR , m 1 " r— (. E 

Xq = hm — : — = lim villi \u\ — — 

R,m^o let dm R,m^o Act \ a 
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Using eqs. (12) and (13), the above equation can be written in explicit form as 




x 



(18) 



This equation depends only on the initial data g and E, as in the dust case, and on the 
choice of the material we are dealing with, i.e. the equation of state h. This means that the 
dynamics has been completely "gauged away" ; such a simplification cannot be achieved in 
comoving coordinates since in such coordinates the general exact solution is not available 
in explicit form. In the next two sections, we shall illustrate a simple way to extract 
physically interesting information from this equation without solving it explicitly. 

3.3 Remarks on shell-crossing singularities 

As recalled above, shell crossing singularities correspond to zeroes of Y', so that in 
the mass-area description a shell crossing occurs when Y jm vanishes. Generally speaking, 
we do not expect a zero of Y^ m to occur at a turning point (u = 0) so that equation (15) 
implies that shell crossing singularities correspond to zeroes of H. Let R sc (m) be the 
curve on which such singularities eventually occur. Using eqs. (13) and (16) we obtain 
that R sc (m) must satisfy to 



This equation may have physically meaningful solutions. For instance, consider the Tol- 
man — Bondi case. The solutions of equation (19) are physically meaningful only if the 
vanishing of Y iTn happens before (in comoving time terms) the singularity at Y = is 
reached. It is possible to characterize fully in terms of differential inequalities the set of 
initial data such that no shell crossing occur in physically allowed "times" (Hellaby & 
Lake 1985, Newman 1986, Jhingan & Joshi 1998). On the contrary, Ori (1990) used the 
charged dust counterpart of this equation (see appendix) to show that the characteristic 
"bounce in a new universe" process (De La Cruz & Israel 1967) which is typical in such 
solutions always occur after a shell-crossing, thereby casting serious doubts on its physical 
realizability (see also Ori 1991). 

In the general case of non-vanishing tangential stresses, the analysis is also possible 
in full generality and will be presented elsewhere. 

IV. The nature of the central singularity: a conjecture 

As we have seen, the nature of the central singularity depends on the existence of 
solutions of equation (18). A complete study of this equation requires a detailed investiga- 
tion on the behaviour of the equation of state in the limit of approach to the singularity in 




(19) 
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physically valid situations, and goes far beyond the scope of the present paper. However, 
some insights into this problem can be obtained by a careful analysis and comparison with 
the (already well known) results holding for dust. For our considerations it will be suffi- 
cient to consider the case of marginally bound collapse; we are, therefore, going to give 
a simple derivation in mass-area coordinates of the results on the nature of the central 
singularity in this case (we completely refer the reader to the original paper by Singh & 
Joshi (1996) for details). 

The marginally bound dust case corresponds to E = h = 1. From formula (16) we get 

g(m) =±J—m*-£-(m-*R°). (20) 

The behaviour of this function as m tends to zero can be obtained as follows. Consider 
regular initial data in comoving coordinates. Then the function F will be of the form 

F(r) = F r 3 + F q r q+3 + ... 

where q is the order of the first non vanishing derivative of the initial density profile at the 
centre, and dots stand for higher order terms. Therefore, we have 

R\m) = F-\m) = ( - A ( ™) ^ + ... 
1 ' ' \F J 3F \F J 

Considering now equation (20) and recalling that we are considering collapse (so that the 
negative sign must be chosen) we obtain 

q_ -t 

g P q m s , 

where 

P ._ _qf* 



9V2F ' 



-+ 2 

3 ' 2 



Thus g exhibits a "critical" behaviour: it diverges (respectively, goes to a finite non-zero 
limit, vanishes) if q < 3, q = 3, q > 3. Surprisingly enough, it is this behaviour that 
governs the nature of the singularity. In fact, equation (18) yields 



%0 = 7T~ lim 



o 



l-a 



' I 



The first term in square brackets goes to a finite, non zero limit iff 

a = - ( 1 + 

so we get 



3 V 3 



1 / 2 A i /mH 1 -!) 

x ° = |(TTf) ( Pg - v^o R ^o i-Tx- " 1 1 • ^ 21 ) 
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If q is "super-critical" (q > 3) the limit diverges: there are no null geodesies escaping and 
therefore the singularity is not naked. If q is "sub-critical" (q = 1,2) the limit goes to 
minus one and (21) gives a real positive solution for xq (provided that the initial density 
is decreasing outwards): the singularity is naked. At the critical value q = 3 (21) becomes 
a quartic equation. This equation has no real positive roots (and therefore the singularity 
is covered) if the quantity ( = F 3 /(2a/2Fq^ 2 ) is greater than a certain numerical value, 
otherwise the singularity is naked. 

The qualitative features of the general (i.e. non-marginally bound) case are similar 
to those recalled above, namely nakedness depends on the "critical" behaviour of some 
parameter q (which reduces to the parameter q in the marginally bound case); if the 
singularity is censored the limit diverges for any q greater than the critical value (Singh & 
Joshi 1996). 

Consider now a generic solution with tangential stresses. To identify it uniquely, we 
need to chose the equation of state and the initial distribution of density and velocity. This 
means that the space of the free functions can be visualized as follows: to any fixed choice 
of the initial data g{m) and E{m) corresponds a family of solutions Sh- Each member 
of this family corresponds to a different material (a different choice of h(m, R) within 
the physically allowed range) and each family contains one and only one Tolman-Bondi 
solution <Si ("dust limit") corresponding to h = 1. Chosen a family Sh, we can immediately 
infer from the Singh & Joshi work if the dust limit Si corresponds to a naked singularity 
or to a blackhole. 

We conjecture that, if the central singularity of Si is not naked and it is not critical 
(i.e., the limit in (21) is divergent), the central singularity of Sh is also not naked for any 
physically valid choice of h. Roughly speaking, this means that one cannot use a physically 
valid tangential stress to undress a covered dust singularity. 

The above conjecture is based on arguments of physical plausibility as follows. Con- 
sider a small deviation from the dust equation of state. This can be represented as 

h = l + n(R,m) , (22) 

where the function fx is positive and vanishes at R = R (m) (a reasonable choice for [i 
could be the "quasi-hookean" equation of state n = /j,o(m)(R — R ) 2 with positive no). 
Then each term in equation (18) can be expanded to first order in \x\ if this function is 
physically valid (i.e. is chosen as described above) every such terms will be quadratic in 
R — R . In particular, the last factor in round brackets will have this behaviour. Now, it 
is easy to check that quadratic terms in R — R° (or higher order terms) cannot regularize 
a diverging behaviour of the zero-order term in this factor. Since divergency occurs in the 
non critical covered case, the conjecture is proved at least for small deviations from the 
dust equation of state. 

V. Self— similar spacetimes 

The above described proposal on the nature of the central singularity, although cov- 
ering a large sector of the space of initial data, leaves completely open the problem of 
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interpretation of critical behaviour. For instance: what happens to marginally bound so- 
lutions having a naked dust limit with q = 3? Do such solutions remain naked with the 
addition of any (physically valid) tangential stress? Do they become always covered? It 
seems likely that none of the above would hold, but rather that the threshold of black hole 
formation for fixed initial data should depend on the equation of state (i.e. on the choice 
of the function h), hopefully in a physically reasonable and understandable way. We don't 
have the answer to this question yet. However, we are going to present here some (again, 
qualitative) evidence that really a behaviour like this should occur. 

Difficulties in studying critical cases arise because we must investigate finite values 
of the limit (18) and, therefore, existence of positive solutions of the root equation. To 
get some insight into this we consider, among the solutions presented above, a particulary 
simple case which, however, is not deserved of physics, namely the case of self-similar 
spacetimes (see e.g. Carr 1997). 

In mass area coordinates we can use as self similar variable the quantity x = R/2m. 
It is easy to check that the spacetimes are self-similar if the following conditions hold (see 
Magli 1993 for a discussion of self-similarity in the case of anisotropic matter): 

1) h is a function of x only 

2) g(m) and E(m) are constant. 

In what follows, it will again be sufficient to consider the marginally bound case E = 1. 
The above conditions imply that R°(m) is a linear function, and indeed it is well known 
that for self-similar spacetimes the mass function (the inverse of R°) is linear. We therefore 
set F(r) = Xr/2 (A = const.) so that 

„ n . . 2m 
R°(m) = — , 

this implies that the value x of x at initial data is 1/A. Using equation (16) with h(x) = 1, 
the root equation (18) can be written as 



1 

X = ~2 



i r 

— + / g(z)dz 

where 7 := A 2 /12 and 




h ( , x dh \ ( 1 1 



3 

5W: =il, 1 + 2 F^Jl,- 1 + F + iJ ' ' (24 » 

To extract from equation (23) some qualitative information we again start from the 
dust limit. What happens in this limit is already well known (Joshi & Singh 1995), and 
we completely refer the reader to this paper for details. 

Setting h = 1 equation (23) simplifies to 



— i + ? J,! h h -1 ■ (25) 
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To facilitate the comparison with the Joshi & Singh paper we change variable to 
then (25) becomes a quartic equation in the variable y: 



y 3 (l-y)= 7(2 - yf , (26) 

where < y < 2/3. Rewriting this equation in canonical form as ay 4 + 4by 3 + 6cy 2 + Ady + 
e = 0, it can be shown that real, positive solutions exist only if the quantity 

A := (ae - 4bd + 3c 2 ) 3 - 27(ace + 2bcd - ad 2 - eb 2 - c 3 ) 2 , 

is negative (in this case there are two such solutions). This happens if 7 < 71 ~ 6.41 x 10 -3 
or 7 > 72 ~ 17.32. Thus the collapse leads to black hole formation if 71 < 7 < 72, to 
naked singularities otherwise. The range 7 > 72 is, however, unphysical since it would 
correspond to imaginary values of x. Therefore self-similar dust spacetimes exhibit a 
"phase transition" between naked singularities and black holes. The transition depends 
on the value of 7 which is the remaining free parameter. This quantity is related to the 
central density of the material; using a fiducial model Joshi & Singh have shown that 7 
typically do belong to the range of black hole formation for densities near the nuclear one. 

We now want to investigate, at least qualitatively, the changes introduced in the above 
picture by the presence of tangential stresses. For, consider once again a small deviation 
from the dust equation of state of the form (22) (obviously, for self-similar spacetimes \i 
has to be be considered a function of x only). In what follows, we shall systematically 
discard terms of order higher than one in fi. Expanding the root equation (23) we obtain 

y 3 (y-f) =i(y-2 + k) 3 . (27) 



3, 

In the above formula, K is the value at order one of the following function 



9 (2 

h {y) U " y j 1 v " 1 ' y i( i- itf ) v ~ V" 1 ' ~~ dz, 

(28) 

Notice that this function is strictly positive. Linearity in fx also implies K = K(yo) where 
yo is the solution of the dust quartic (26) in the neighbourhood of which we want to study 
the deviation from the dust case. 

Expanding also equation (27) to first order, we finally obtain a quartic with "dis- 
placed" parameters a, 6, c, d, e. A quite long, but straightforward, calculation gives that the 
range of black hole formation 71 < 7(< 72) is altered by the perturbation as 71 < 7(< 72), 
where 

71 = 71 - Kh , (29) 

72 = 72 + K82 
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and Si ~ 17 x 10~ 3 , S 2 ^ 12.8 (as usual in any perturbative approach, the above results 
also give the condition of applicability of the approximation: K must be not greater than 
about 1/3). 

Formulae (29) show, at least at a qualitative level, that black hole formation is facili- 
tated by the presence of tangential stresses. Indeed, since K is positive, the upper bound 
72 becomes higher and certainly remains unphysical, while the lower bound 71 tends to 
decrease. It might happen that the addition of tangential stress dresses the singularity, 
which would be equivalent to 71 < 0, but of course to draw a conclusion of this kind it will 
be necessary to investigate the root equation without approximations. 

VI. Concluding remarks 

The results of the present paper can be summarized as follows. First of all, we have 
shown that Ori's mass-area formalism can be used to bring the general spherically sym- 
metric solution of the Einstein field equations with non-vanishing tangential stresses in a 
very explicit form, in which only the calculation of an integral remains to be performed 
(this can be done independently whether electromagnetic coupling is present or not; the 
generalization to charged materials is indeed straightforward and is briefly reported in the 
appendix) . 

The introduction of mass-area coordinates proves to be a very powerful tool as far 
as the analysis on existence and nature of singularities in such solutions is concerned. In 
fact it allows to obtain the root equation governing the nature of the central singularity in 
explicit form. 

We presented first results coming from the investigation of the dust limit of this 
equation. Such results give some insight about what should be the nature of the final fate 
of collapse with tangential stresses. In particular, we proposed a conjecture which, roughly 
speaking, asserts that "tangential stress cannot undress covered dust" . In the last section, 
we discussed self-similar spacetimes and showed, at least qualitatively, that the effect of 
the stress can be an enlargement of the blackhole initial data space. 

Both the above recalled results depend on the structure of the state equation, and 
therefore show once again and in a clear way that a connection should exist between a (still 
lacking) mathematically rigorous formulation of cosmic censorship and the conditions of 
physical acceptability of the equations of state. Such conditions obviously include the 
energy conditions but also the existence of an absolute minimum of the internal energy, 
which is intimately related to stability issues; a relevant improvement in our understanding 
of this topic could come from the knowledge of the explicit structure of the blackhole 
threshold in terms of the derivatives of the state equation evaluated near the singularity. 
Work in this direction is now in progress. 
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Consider a material carrying a no n- vanishing charge density a. We keep the descrip- 
tion of the mechanical and gravitational degrees of freedom as in the body of the paper, 
and simply introduce the Maxwell tensor which, due to spherical symmetry, has only one 
independent component F mR := S(m,R). Maxwell's equations yield 

Q 



£ = 



4ttR 2 VHu 

where Q = Q(m) is arbitrary. The field equation (8) now has a source term and reads 



r 2 Vh ' 

Q ,m 



1 

R 2 



A 



1- — -R — 



H 



A 



H 



.R 



R* ' 



while (9) and (10) remain unchanged. Since Q = Q(m), the equation above can be 
integrated at once and gives 



R R 2 

Using B 2 = H + A/u 2 and equation (7), from (9) we get 



u = ±1 



2m Q 2 1 



R 



R 2 h 2 



Q 



R 



where the "specific charge" £ := EQ :Tn . The metric function B can be written as 



hu 



R 



H 



so that (10) gives 



It follows 



h 



y/H(m, R) = g(m) ± J G(m, R)dR , 
where g(m) is arbitrary and 



G := 
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For h = 1 (dust case) the function G is a rational fraction and its integral may be carried 
out explicitly (Ori 1990), while for Q = (non-charged case) the above function coincides 
with that defined in (14). 



